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The specific heat of superfiuid ^He, disordered by a silica aerogel, is found to have a sharp 
discontinuity marking the thermodynamic transition to superfluidity at a temperature reduced from 
that of bulk ^He. The magnitude of the discontinuity is also suppressed. This disorder effect can 
be understood from the Ginzburg-Landau theory which takes into account elastic quasiparticle 
scattering suppressing both the transition temperature and the amplitude of the order parameter. 
We infer that the limiting temperature dependence of the specific heat is linear at low temperatures 
in the disordered superfiuid state, consistent with predictions of gapless excitations everywhere on 
the Fermi surface. 

PACS numbers: 67.57.-z, 67.57.Bc, 67.57.Pq 



Two essential characteristics of superconductors are 
manifest in the specific heat. First, there is a jump at the 
transition temperature of magnitude that depends on the 
strength of the electron pair interactions that lead to su- 
perconductivity. Secondly, the temperature dependence 
at low temperature gives a fingerprint of the energy gap 
structure. Consequently, the measurement of the specific 
heat is one of the most fundamental to understanding the 
nature of superconductivity. In fully gapped supercon- 
ductors, such as aluminum, the specific heat decreases 
exponentially with decreasing temperature following an 
Arrhenius relation Ij, thereby demonstrating unambigu- 
ously the existence of a full gap in the electron quasi- 
particle excitation spectrum and providing a direct mea- 
surement of its size. Both of these basic thermodynamic 
behaviors conform to predictions of the Bardeen, Cooper, 
and Schrieffer theory 2] and can be expected to hold in 
general for any condensate of fermion pairs into a state 
that is fully gapped. This is precisely what is expected 3] 
and observed 4J for the B-phase of superfiuid ^He, a p- 
wave superfiuid with an isotropic energy gap, similar to 
conventional 5- wave superconductors. 

Shortly after the success of BCS theory it was sought 
to understand how perturbations, such as impurities, 
modify superconducting behavior. It was found that 
magnetic scattering of quasiparticles suppresses both the 
transition temperature and the gap magnitude leading 
to 'gapless superconductivity' 2]. The same should be 
true for all superconductors and fermion superfluids even 
those that do not have 5-wave states in which case all 
forms of elastic scattering will produce these suppression 
effects p. In particular it should be true for superfiuid 
^He. In this letter we present systematic measurements 
of the heat capacity in the B-phase of superfiuid ^He, dis- 
ordered by impurities, demonstrating suppression of the 
transition temperature, suppression of the order param- 
eter, and gapless superfiuid behavior. 

Superfiuid ^He is a unique example of a Cooper pair 
condensation. It was the first unconventional pairing 



state discovered. By unconventional we mean that there 
are symmetries of the normal Fermi liquid that are spon- 
taneously broken in the superfiuid phases in addition to 
gauge symmetry, e.g. spin and spatial rotations. In this 
sense ^He is similar to cuprate superconductors and to 
the heavy fermion superconductor UPta as well as a num- 
ber of other recently discovered superconducting com- 
pounds. However, in contrast to the superconductors, 
superfiuid ^ He is a neutral system, easily prepared as the 
purest known material, whose general properties are very 
well established. It is not so straightforward, however, to 
introduce impurity scattering in a controlled way in ^He. 

Porto et al. and Sprague et al. j^, |2| found that they 
could use very dilute silica aerogels to introduce impurity 
scattering in superfiuid ^He and subsequently there has 
been a wealth of new information about the nature of this 
modified superfiuid state. Impurity scatte ring models for 
^He were developed by Thuneberg et a/.[l^ with clear 
predictions for the suppression of the transition temper- 
ature and the heat capacity discontinuity. Sharma and 
Sauls ^3 have calculated transport behavior as well as 
the density of quasiparticle excitations inside the gap 
which define its gapless superfiuid behavior. Our ex- 
perimental results are quantitatively consistent with the 
calculations of these basic thermodynamic properties. 

Earlier work of He et a/.p^] and preliminary results 
from our laboratory have provided convincing evi- 
dence for the existence of a heat capacity anomaly at the 
transition temperature. Additionally, Fisher et a/.fl^ 
found that the thermal conductivity in the disordered 
A-superfluid state in a magnetic field has a linear tem- 
perature dependence when extrapolated to very low tem- 
pertature. If the heat current were carried by gapless ex- 
citations at low energy there should be evidence for these 
excitations in the low temperature specific heat. 

Our low temperature calorimeter was cooled with a 
PrNis adiabatic nuclear demagnetization refrigerator to 
temperatures below 0.75 mK above which we made mea- 
surements of heat capacity using an adiabatic calorimet- 
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ric method. Temperature was determined from mea- 
surement of the ac susceptibiUty of a lanthanum diluted 
cerium magnesium nitrate salt, calibrated with respect to 
a ^He melting curve thermometer according to the Grey- 
wall temperature scale*?]. The calorimeter volume was 
measured to be 1.91 ± 0.02 cm^ with the aerogel sam- 
ple in place, leaving 44 % of the space available for bulk 
^He. The porous silica aerogel was in the form of a right 
circular disk 0.380 cm thick and 1.902 cm in diameter. 
From its dry weight, 42.4 mg, we found the porosity to 
be 98.2%. The sample was first cooled by nuclear de- 
magnetization. By opening a cadmium superconducting 
heat switch we could isolate the sample cell from the nu- 
clear stage, allowing the ^He to warm up slowly owing 
to the ambient heat leak, typically Q ~ 0.1 nW, starting 
in the B-phase 19j. Heat pulses were applied with AQ in 
the range 0.1 ~ 0.5 /iJ for 30 ~ 60 s, and the temper- 
ature increments AT were measured to obtain the heat 
capacity Ctot = AQ/AT shown in Fig. 2 after subtrac- 
tion of the addendum. A warm-up trace without heat 
pulses, in the form = C/Q, is shown in the inset to 
Fig.l providing our highest resolution of both transition 
temperatures, in bulk liquid ^He and ^He in aerogel. 

We assume that liquid ^He above its transition to su- 
perfluidity is a Fermi liquid with known interaction pa- 
rameters even within the aerogel. This is reasonable since 
the aerogel structure is on a much larger scale than the 
Fermi wavelength 20j. Using the specific heat measure- 
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FIG. 1: Pressure temperature phase diagram for bulk ^He 
and disordered ^He in aerogel. The smooth curve is an 
interpolation from Greywall's data0 for bulk ^He. The 
data points are high resolution measurements of the tran- 
sition temperatures obtained during slow warming shown 
in the inset at P = 20.06 bar. Two heat capacity anoma- 
lies are evident. The solid fit curve defines the bulk tran- 
sition (dotted vertical line, open circles). The superfluid 
transition temperature in aerogel is at a lower tempera- 
ture (solid vertical line, solid circles), ^ 80/iK wide. 
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FIG. 2: Heat capacity of ^He in aerogel. On the left side 
we show the heat capacity as a function of temperature 
for bulk ^He and ^He in a 98% porous silica aerogel con- 
tained together in the same calorimeter with addendum 
subtracted. On the right side we show the heat capacity 
divided by the temperature for the disordered superfluid 
alone. The curves represent the known temperature de- 
pendence of the heat capacity for bulk ^He (dashed) and 
a smooth fit to the data in the disordered superfluid with 
the constraint for entropy conservation (dotted) described 
in the text. 



ments by Grey wall on bulk ^He we can deduce from 
our experiments, a) the calorimeter background at tem- 
peratures above the superfluid transition in aerogel, b) 
the volumes of ^He in the aerogel and in the bulk, c) the 
heat capacity jump in the disordered superfluid phase, 
and d) the temperature dependence of the heat capacity 
in the disordered superfluid state. 

The background heat capacity, Cadd^ can be entirely 
attributed to solid ^He on the surface of the aerogel. Two 
layers of ^He are paramagnetic jol [Tsl] and can be removed 
by substitution with two layers of non- magnetic ^He re- 
ducing the addendum to less than 0.25 mJ/K. The helium 
isotope mixture experiments were limited to > 4 mK ow- 
ing to long thermal time constants in the calorimeter. 
The temperature and pressure dependence of the adden- 
dum is shown in Fig. 3, of a comparable magnitude 
to that reported in vycor glass [l^. 

The measured heat capacity, Ctot^ corrected for the 
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addendum, Cadd, is 



C(T, P) = Ctot - Cadd = VaCa + VbCb 



(1) 



shown in Fig. 2. Above the bulk transition tempera- 
ture, Tc, C{P,T) is proportional to temperature, since 
the specific heats of the fluid inside, c^, and outside, 
c^, the aerogel are identical: Ca = cb = where 
7 = {tt'^ /3)k%No and Nq is the density of states at the 
Fermi energy. Knowing the total volume Va + Vb of aero- 
gel and bulk fluids, we can determine the addendum in 
this temperature range. At Tc the bulk heat capacity 
jump was compared with Grey wall's values to give a di- 
rect measure of Va from which we found the fluid volume 
in the aerogel to be 1.028 ± 0.021cm^, consistently at all 
pressures and close to the measured geometry of the 
aerogel giving a pore volume of 1.062 ± 0.006 cm^. With 
Va and Eq. 1, we can determine the addendum in the 
range Tc > T > Tc^a^ by comparing our measurements 
with the known specific heat of the bulk superfluidy. 
We must also allow for contributions to the heat capac- 
ity from a normal state region at the surface of the heat 
exchanger, approximately a coherence length, ^(T), in 
depth. This small contribution is found by fitting the 
slow warm-up curves as shown in the inset to Fig.l to ob- 
tain an effective area of the heat exchanger A = 1.71m^, 
valid at all pressures. The temperature dependence of 
the addendum at 1.02 bar is smooth to below 1 mK. 
Consequently, it is reasonable to assume that this is the 
case at higher pressures and to extend linearly the tem- 
perature dependence of the addendum from Tc,a to lower 
temperatures in order to extract the heat capacity, Ca, 
of the disordered superfluid ^He, Fig. 2. 

We show Ca/T on the right side of Fig. 2 after subtract- 
ing the bulk contribution, the surface contribution, and 
the addendum. The heat capacity jumps appear quite 
sharp in the figure but substantially reduced even com- 
pared to the weak coupling BCS value of AC/C = 1.43. 
They are shown as a function of pressure in Fig. 4, nor- 
malized to their bulk values. Consequently, the ampli- 
tude of the order parameter is also reduced since the 
jump is related to the temperature dependence of the 
gap function, Aa(T), near T^a, in the Ginzburg-Landau 
limit. 



Aa{T) = 7rkBTc,a[ 
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At lower temperatures, below Tc,a, we note that Ca/T 
appears to be quite linear. We take advantage of this 
fact and we invoke the third law of thermodynamics. 
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to determine the limiting behavior of as T ^ shown 
by dotted curves in Fig. 2. Graphically, this means that 
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FIG. 3: Heat capacity addendum. The temperature and 
pressure dependence of the heat capacity addendum is 
shown for bulk ^He in aerogel with an estimated accuracy 
of 0.25 mJ/K. 



for each pressure the two shaded areas must be equal. If 
the heat capacity is a monotonic function of temperature, 
as is expected theoretically, then to satisfy Eq. (3), Ca/T 
must have a non-zero intercept at T = requiring a non- 
zero density of states, A^a(O), at the Fermi energy in the 
superfluid state. The curves describing this behavior in 
the disordered superfluid were constructed phenomeno- 
logically to satisfy the constraints and our interpretation 
does not depend on its functional form nor is it particu- 
larly sensitive to the data at the lowest temperatures. For 
example in Fig. 2 at P = 29.02 bar the data fall below the 
curve likely owing to our overestimate of the addendum. 
The low energy excitations affecting the low temperature 
specific heat are a broad band of Andreev bound states, 
centered at the Fermi energy, that form near the aerogel 
strands llj. As a result superfluid ^He-B in aerogel is a 
'gapless superfluid'. 

The simplest theoretical description 0, considers 
isotropic quasiparticle scattering spread homogeneously 
throughout the aerogel volume, characterized by a sin- 
gle parameter, the transport-mean- free path. At- This 
model (HISM) provides a qualitatively consistent account 
of the superfluid transition temperature suppression [loj. 
the suppression of the A-phase and the po lycr itical 
point 'l^ measurements of spin diffusion and 

the thermal conductivity 111 Here we see that the 
model gives a quantitatively consistent picture for the 
heat capacity jump over a wide range of pressure. A 
calculation lOj with the HISM gives a best fit value 
of At = 180 nm. In this fit we rescale strong cou- 
pling contributions jlSl to the free energy by the ratio 
Tc,a/Tc where Tc,a is calculated self-consistently within 
the HISM. We also find reasonable agreement with the 
prediction from the HISM for the density of states of 
gapless excitations. The calculations as a function of 
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FIG. 4: Heat capacity jump and density of states. The 
pressure dependence of the heat capacity jump (open cir- 
cles, left axis) and the density of states (filled circles, 
right axis) at zero energy (the Fermi energy) is shown 
for bulk ^He in a 98% porous silica aerogel. The error 
bars arise from possible systematic error in the choice of 
temperature dependence for the heat capacity just below 
Tc^a- Other statistical errors are of the size of the data 
points. Calculations from the HISM are shown as curves 
for various values of the transport mean free path. 

pressure for A/'a(0) are shown in Fig. 4 for several values 
of At near 150 nm and are qualitatively consistent with 
our analysis of the heat capacity jump, At = 180 nm. 

Discrepancy between the experiment and the model 
for the density of states can be attributed in part to es- 
timated systematic errors including the extrapolation of 
the addendum to T < Tc,a and possible anisotropic or 
inhomogeneous scattering. Nonetheless, the consistency 
that is apparent within the heat capacity experiment pro- 
vides strong evidence that elastic scattering suppresses 
superfluidity and is approximately isotropic and homo- 
geneous. It is interesting to compare our results with 
unconventional superconductors in several cases. Preci- 
sion measurements of the transition temperature of UPta 
indicate suppression of the transition temperature 
with increased scattering but, in this case, with a sig- 
niflcant anisotropy. In cuprates doped with Zn impurity, 
such as YBa2(Cui_^Zn^)307, suppression effects have 
been observed in the heat capacity strikingly simi- 
lar to those we report here. 
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